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We propose a technique for robust and efficient naviga- 
tion in the Hilbert space of entangled symmetric states of a 
multiparticle system with externally controllable linear and 
nonlinear collective interactions. A linearly changing exter- 
nal field applied along the quantization axis creates a network 
of well separated level crossings in the energy diagram of the 
collective states. One or more transverse pulsed fields ap- 
plied at the times of specific level crossings induce adiabatic 
passage between these states. By choosing the timing of the 
pulsed field appropriately, one can transfer an initial product 
state of all N spins into (i) any symmetric state with n spin 
excitations and (ii) the JV-particle analog of the Greenberger- 
Horne-Zeilinger state. This technique, unlike techniques using 
pulses of specific area, does not require precise knowledge of 
the number of particles and is robust against variations in the 
interaction parameters. We discuss potential applications in 
two-component Bose condensates and ion-trap systems. 



I. INTRODUCTION 

Entanglement is a unique quantum feature which has 
enjoyed considerable attention in the last few years. It 
plays a crucial role in many rapidly developing areas of 
contemporary quantum physics, such as quantum infor- 
mation JIJ and fundamental tests of quantum mechanics 
Various quantum systems have been suggested for 
controlled creation of entanglement, e.g. trapped ions 
||, spins in magnetic field Q, quantum dots ||, cavity- 
quantum-electrodynamics systems || , crystal lattices , 
Josephson junctions ||, and others. 

In order to entangle N spin-i particles interaction be- 
tween the spins is required and external control of this 
interaction is necessary to generate specific many-particle 
states. For the latter purpose one can use sequences of 
resonant external pulses of precise area, e.g. 7r-pulses. 
While this technique is conceptually simple, it is very 
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sensitive to variations in the pulse area and resonance 
mismatch, which can be caused by temporal and spatial 
fluctuations of the external field and may lead to signif- 
icant errors. Hence an important practical challenge is 
to design robust and efficient methods for a controlled 
navigation in the multi-particle Hilbert space. This is of 
particular importance in mesoscopic systems, where only 
limited control over the interaction parameters and the 
number of particles is possible. 

The simplest interaction that can lead to entanglement 
in a collection of spins involves either pairwise nearest- 
neighbor interactions or a collective coupling between all 
particles. An example for the first case is the Ising model 
Q , while the collective coupling of ions to a phonon mode 
in an ion trap || or the self-interaction in a Bose- Einstein 
condensate (BEC) jlO| is an example for the second. 

We here analyze the second type of systems and pro- 
pose a method for controlled, efficient and robust naviga- 
tion in the space of symmetric multi-particle entangled 
states. No precise knowledge of interaction parameters or 
particle number is needed and only certain adiabaticity 
criteria have to be fulfilled. The proposed technique is a 
multi-particle generalization of our earlier proposal EJ] 
for creation of entanglement in a pair of two-state sys- 
tems by using adiabatic passage induced by a suitably 
crafted external field. The multi-particle problem adds 
some new challenges as it involves in general multi-step 
as well as direct transitions; this opens a variety of paths 
between any pair of multi-particle states. We use this 
to advantage and demonstrate that certain paths are in- 
sensitive to the number of particles N , a property that 
is particularly significant for systems with large N, for 
which N is not usually known exactly. 

The scheme proposed in the present paper allows to 
create robustly and efficiently maximum entanglement 
starting from a product state. In particular, it can be 
used to create various special entangled states, such as of 
the Greenberger-Horne-Zeilinger (GHZ) type or the so- 
called IF-states. Moreover arbitrary transitions among 
the collective states can be realized. After introducing 
the concept of the method we discuss two specific imple- 
mentations: an ion-trap scheme similar to that of M0lmer 
and S0renson jl2| and a coupled two-component BEC in 
the two-mode approximation paj. 
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II. NAVIGATION IN THE HILBERT SPACE OF 
THE COLLECTIVE STATES 

A. Controlled collective spin interaction 

We consider a collection of N identical spin-i particles 
with a total angular momentum operator J. The simplest 
collective interaction that allows to entangle individual 
spins is quadratic in one of the components of J, viz. (in 
units h = 1) 



H(t)=U 2 z +B(t)-J, 



(1) 



where £ is the spin-spin interaction constant and B(f) 
is some time-dependent external field. Special cases of 
this type of Hamiltonian have been discussed by several 
authors (see, e.g., [O and references therein). 

Since the Hamiltonian (|I]) contains only collective spin 
operators, only multi-particle states with the same sym- 
metry upon particle exchange are coupled. In particu- 
lar, if the system is initially in the product state of all 
spins in the spin-down state |44 ... J.), only symmetric 
collective states with the same (maximum) angular mo- 
mentum J = jN will interact. Thus the 2 Ar -dimensional 
Hilbert space reduces to the (N + l)-dimensional sub- 
space of symmetric multiparticle states. Because each of 
these states is characterized by a definite number of exci- 
tations n = 0, 1, . . . , N (which corresponds to an angular 
momentum projection m = — J + n G { — J, J}), we shall 
use n to label the states. They are given by 
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where af is the Pauli spin-flip operator, which inverts 
the spin of the ith particle. In particular, the single- 
excitation symmetric state |1) is an iV-particle analog of 
the W -state p|: 

|i) = -4={ ITU ■ • • I) + 1444 ■ • • 4H — llll •••!)}• 



N 

This state is maximally robust against disposal of any of 
its qubits Obviously, state \N — 1), which has N — 1 
spins up and one spin down, is also of the VF-type. 

We assume, similarly to , that the particle-field cou- 
pling consists of a linearly changing component along the 
z-direction (chosen as the quantization axis) and a pulsed 



component in the x-direction. Then the Hamiltonian ([[]) 
takes the form 



H(t)=£J* -AtJ z +Q(t)J x , 



(3) 



where A is assumed real and positive. Because the opera- 
tor J x connects only adjacent states, the linkage pattern 
is chainwise. The direct coupling between each pair of 
adjacent states \n) and \n + 1) is given by 



n 



n,n+l 



(t) = ^J(J +1) -m(m + l)Q(t), 



(4) 



i.e. all couplings have the time dependence of the exter- 
nal pulse, but different magnitudes. 

The energy of each state \n) changes linearly in time 
with a slope proportional tom = -J + n, 



E n (t) = rn 2 £ - mAt. 



(5) 



This creates a web of level crossings in the energy dia- 
gram of the collective states as shown in Fig. [j] for the 
case of N = 4 particles. These level crossings can be used 
to design various navigation routes between the collective 
states. 

As follows from Eq. (||), states \n) and |fe) cross at 
time 



tnk = {n + k-N)j. 



(6) 



These crossings are equidistant and separated by a time 
interval r = £/A, which does not depend on the number 
of particles N, but only on the interaction parameters. 
Pairs of states with the same total number of excita- 
tions n + k cross at the same time. However, due to the 
presence of the nonlinear interaction term m 2 £, all level 
crossings are well separated in energy. 

The positions of the crossings may or may not depend 
on N. For example, the crossing between states |0) and 
|1) is situated at ioi = — (N — l)£/A. Hence this time 
will be known only if N is known exactly; we will see be- 
low that this restricts the possible scenarios of adiabatic 
transfer. In contrast, the crossing between the two prod- 
uct states |0) and \N) does not depend on iV, because it is 
situated at iow — 0, and the time t = is determined by 
the zero value of the longitudinal field, B z (0) = 0. The 
crossing between states \N) and |1), which is situated at 
tiN — S./A, does not depend on N either. 



B. Principles of navigation in Hilbert space 

1. Navigation routes 

Once a network of level crossings is created, one can 
design in principle any navigation route in the Hilbert 
space by choosing properly the timing of the pulsed ex- 
ternal field. It is most natural to assume that the multi- 
particle system is prepared initially in one of the product 
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FIG. 1. Energies of the five eigenstates of the Hamilto- 
nian (|^) for TV = 4 particles. At equidistant times t n k there 
are diabatic level crossings between the energies (^) of the 
collective states. States |0) and |4) are product states, all 
others are entangled symmetric states. Application of an ex- 
ternal coupling pulse at time t n k leads to an avoided crossing 
and adiabatic population transfer between states |n) and \k). 
The four frames apply to the cases when an external coupling 
pulse is applied at the crossings (a) toi, (b) to2, (c) to3, (d) 
toA- As a result, four navigation routes are created connecting 
the initial state |0) to states |1), |2), |3), and |4), respectively. 




FIG. 2. Numerically calculated populations of the sym- 
metric collective states \n) in a four-particle system, initially 
in the product state |0), plotted against the center To of a 
Gaussian coupling pulse, Q(t) = S7oexp[— (t — Tb) 2 /T 2 ], for 
£ = 20T~\ Q = 50T -1 and A = 5T~ 2 . 



the range of times over which complete transfer occurs 
decreases as the number of excitations associated with 
the transfer increases. 

The present technique allows also to transfer entan- 
glement. Indeed, if the many-particle system is initially 
in the entangled state |n), it can be transferred adia- 
batically to another entangled state |fc) by applying an 
adiabatic pulse at the crossing time t n k of these states. 
An entangled state \n) can also be transferred into one 
of the unentangled states |0) or \N); this can be used for 
measurement of entanglement. 



states, for example in |0). If the system evolves along this 
state and the designed route requires that it must make a 
transition to the entangled state \n), one should apply a 
sufficiently strong (adiabatic) external pulse at the time 
of diabatic crossing to n between |0) and \n). The interac- 
tion will open an avoided crossing between the energies of 
the corresponding adiabatic states and will force the sys- 
tem to make an adiabatic transition from |0) to \n) [ |L6| . 
On the contrary, if the route requires the system to re- 
main in the same collective state, one should ensure that 
there is negligible interaction between this state and the 
other collective states. With the leeway in choosing the 
time dependence and the intensity of the pulsed external 
field, one can link any initial state to any final state by 
using one or more suitably timed pulses. 

In Fig. ^ we show the numerically calculated popu- 
lations of the five symmetric collective states \n) in a 
four-particle system plotted against the center To of a 
Gaussian coupling pulse. Depending on the timing of 
this pulse, the population is transferred from the initial 
state |0) to different collective states. In agreement with 
the above discussion, the maximum transfer efficiency for 
each collective state \n) is achieved when To is near the 
crossing to n between |0) and |n). One also recognizes that 



2. Conditions 

For such a state engineering to be successful, the level 
crossings must be well separated, i.e. the width T of the 
external pulsed field must be small in comparison with 
the time separation between the crossings: T <C £/A 

On the other hand, the condition for Landau-Zener 
population transfer around the crossing at t n k can be 
shown to lead to the following conditions 

n nfe (i„ fc )T» VlT»l, (7) 

where tt n k(t n k) is the effective coupling between states 
\n) and |fc), estimated at the crossing time t n k- For ad- 
jacent states, which are connected directly, this coupling 
is given by Eq. (Q). For example, the above conditions 
suffice to estimate the feasibility of the direct transition 
10) - |1>- 

For states that are coupled via one or more interme- 
diate states, Q n k(t) is an effective multi-quantum cou- 
pling. This coupling can be estimated perturbatively 
when <C N£ [jlTj by eliminating adiabatically the off- 
resonant intermediate states, which yields 

/ n xl«-fc| 
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FIG. 3. Alternative routes for the transition between the 
product states |0) — > |4) in a four-particle system: a single 
narrow pulse applied at the crossing to4 between |0) and |4) 
(top), and a wide pulse covering all crossings (bottom). 

Thus for £1 <C N£ the effective coupling between \n) and 
|fc) is very small and cannot induce adiabatic evolution. 
Thus we have to consider the case fi > i\T£, which is 
however accessible only numerically. 

C. Choice of navigation path 

As the energy diagram in Fig. [l] suggests, there are 
multiple paths linking each pair of collective states. Each 
of these paths has certain advantages and disadvantages, 
depending on the particular experimental situation. 

1. Transition between the product states 

The two product states |0) and |A") can be linked in 
several different ways. First, one can apply a single adi- 
abatic pulse at their level crossing at time toN — 0, as 
shown in Fig. [| (upper frame) . This approach is indepen- 
dent of the number of particles N because the crossing 
time toN — is well defined by the zero of the linearly 
increasing z-field. In other words, even if we do not know 



the exact N, we can find the exact location of the cross- 
ing between the completely unexcited state |0) and the 
completely excited state \N). However, using this cross- 
ing between |0) and 1 7V) requires much stronger field be- 
cause these states are not coupled directly, but only via 
an A^-quanta transition. 

Alternatively, one can use a train of pulses centered 
at each crossing, so that the population will flow from 
|0) through all intermediate states to reach \N) at the 
end: |0) -> |1) ->• |2) -> ... \N). Because this 
navigation route (the lowest solid curve in Fig. |^) 
passes only through crossings of directly coupled adja- 
cent levels, much less field intensity is needed to satisfy 
the adiabatic condition. However, the first crossing at 
ioi = — (N — l)£/A depends on the number of particles 
N, i.e. this approach is only applicable if A" is known 
exactly. 

A third possibility, applicable only to the transition 
|0) — > \N), is to apply a sufficiently wide single pulse, 
covering all crossings, as shown in Fig. [| (lower frame). 
In this case, only an approximate knowledge of N is re- 
quired, because there are no stringent restrictions on the 
pulse width and timing. 

2. Creation of W -states 

One faces similar choices for the transitions from one 
of the product states to any entangled state \n). For ex- 
ample, the A^-particle W-state can be created using the 
transition |0) — > |1) by applying a single adiabatic pulse 
at time tax = —(N — 1)£/A, as shown in Fig. ^ (upper 
plot). Because this is a transition between adjacent col- 
lective states, only a moderately strong field is required. 
However, this approach is only applicable if the num- 
ber of particles is known precisely because the crossing 
time ioi depends on A^: an error in N, even as small 
as AA^ = 1, is inadmissible since then the pulse will be 
applied at a wrong crossing. 

Alternatively, state |1) can be populated using the 
transition \N) — » |1), as shown in Fig. [| (lower plot). 
For this path, the crossing time £jyi = £/A does not de- 
pend on N. However, because this is an (N — l)-quanta 
transition, the coupled states differ by a large number 
of spin excitations. Then a much stronger field may be 
needed to widen the much narrower avoided crossing and 
induce adiabatic evolution. 

Similar conclusions apply to the other W-state |Af— 1) 
and to any other entangled state \n). 

3. Feasibility of the N -invariant scenario 

The AMnvariance of the latter approach, which uses 
the (multi-quanta) crossings near t — rather than the 
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FIG. 4. Alternative routes for the creation of the W state 
1) in a four-particle system: |0) — > |1) (upper plot) and 
4) — » |1) (lower plot). 

(single-quanta) outside crossings, is a very attractive fea- 
ture because it allows to use this technique without know- 
ing the precise number of particles N. The only problem 
is that for multi-quanta transitions a stronger external 
field is needed to induce adiabatic passage. Therefore 
we have performed numerical simulations to estimate the 
minimal pulse area *4 m i n needed to achieve 90% transfer 
efficiency for the transition \N) — > |1) from the product 
state \N) to the VF-state |1). This area is plotted in Fig. || 
as a function of the particle number N. 

As the figure demonstrates, A m i n increases nearly 
quadratically for small TV and approaches a linear depen- 
dence for large N. This (slow) linear increase implies that 
the conditions on the required resourses for application 
of this technique to many-particle systems are not very 
strong. The linear behavior can be understood qualita- 
tively by noting that, as the numerical simulations show, 
the coupling term J x in Eq. (|J) (which is cx Slo oc ^4 m in) 
is of the same order of magnitude as the vertical energy 
splitting (which is cx Nt;). Hence, the transition proba- 
bility must scale with the parameter £Iq/(N£). 
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FIG. 5. Numerically calculated minimal pulse area ,4 m m 
(defined as the pulse area for which 90% transfer efficiency is 
achieved) for the transition from the product state \N) to the 
IF-state |1) as a function of the particle number N. We have 
used a Gaussian pulse, Q(t) = flo exp [-{t - T ) 2 /T 2 ], with 
£ = 20T~ 1 ,A = 10T" 2 ,T = £/A 

D. Generation GHZ states 

Another interesting application of the present scheme, 
after a slight modification, is the creation of the so-called 
GHZ state @, 

|GHZ} = ^(|0) + |AO). (8) 

Starting from the product state |0) one has to create 
first an equal superposition of this state and the VF-state 
|1), for example, by a single 7r/2-pulse. Then the present 
technique can be used to transfer the population of 1 1) to 
the other product state \N) by applying a single adiabatic 
pulse at the crossing tijv, 

|0>^^(|0) + |l)) ad ^ b ^(|0) + |iV)). 

III. IMPLEMENTATIONS 

In the following we briefly discuss two implementations 
of the technique discussed above. The first one involves 
two-photon transitions in trapped ions, the second a co- 
herently coupled two-component Bose-Einstein conden- 
sate. 

A. Ion-trap system 

Recently S0renson and M0lmer suggested a realization 
of a nonlinear Hamiltonian in a two-level ion trap sys- 
tem [Q, displayed in Fig. || In their scheme, two laser 
fields (in the Lamb-Dicke limit of light coupling) are ap- 
plied with frequencies u>\ and u>2 that are symmetrically 
detuned from the single-photon resonance by a detuning 
5. Then there is a two-photon resonance between states 
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FIG. 6. Ion-trap system for realization of an effective non- 
linear spin Hamiltonian. The application of a symmetrically 
detuned bichromatic field to ions with ground state \g) and 
excited state |e) leads to two-photon resonance between the 
collective states \gg) and |ee) independent of the trap oscilla- 
tion quantum number n. 

with two ions in the ground level |<7<?) and two ions in the 
excited level |ee). If the detuning is larger than the Rabi 
frequencies of the two lasers (5 ^> f2i = O2 = ft) there is 
no single-photon excitation. 

In a coarse-grained time-averaged picture an effec- 
tive Hamiltonian proportional to £ J 2 emerges [fl2f , with 
£ = 2rj 2 Vl 2 v / {v 2 — (5 2 ), 77 being the Lamb-Dicke param- 
eter, and v the phonon frequency. If, in addition to the 
detuned bichromatic fields, a resonant laser coupling the 
z component and one coupling the x component are ap- 
plied, the Hamiltonian is obtained. Since £ is pro- 
portional to the Lamb-Dicke parameter rj ~ 1/ViV, the 
crossing times (^) depend on the number of particles. 
However, this is not a significant problem because N is 
usually precisely known for trapped ions. 

B. Two-component BEC 

Another example for our model Hamiltonian (^|) is the 
trapped atomic BEC in two different internal states a 
and b |I(|[l3|], demonstrated in recent experiments with 
the hyperfme states \F = l,M = ±1) in sodium |l9| . 
The external coupling field B(t) can be provided by a 
chirped radio-frequency pulse. 

The Hamiltonian of a two-component BEC with s- 
wave scattering and coherent coupling between the com- 
ponents in the two-mode approximation is given by 

H{t)=H + H int (t), (9a) 
H = E a a^a + E b tfb + ^a^ 2 a 2 + ^6+ 2 b 2 , (9b) 

H mt (t) = ^oW& + fi(i)[a6 t e-^( t ) + brfe**®]. (9c) 

Here a and b are the annihilation operators for the bosons 
in the two states \a) and |6), and E a and Eb are their 



energies. U aa , Ubb and U a b characterize the particle- 
particle interaction when the particles are in states \a) 
and |6). The time-dependent functions f2(i) and <f)(t) 
are the amplitude and the phase of the coherent cou- 
pling between the two components. If the elastic s- 
wave scattering lengths of the two components are equal 
(U aa = Ubb = U), as in the sodium experiment |l9| |, 
the nonlinear spin coupling £ will be independent of the 
number of particles N. 

In terms of the Schwinger bosonic represenation of the 
angular momentum, 

j x = i(a6 t +6a t ), J y = -\i{atf -6a t ), J z = | (0*0-6+6), 

the Hamiltonian (j|) takes the form 

H(t) = aJ z + £J 2 + Sl(t)[J x cos (f>(t) + j y sin </.(*)], 

where a — E a — Eb and £ = U — \U a b- After a unitary 
transformation ^(t)) = expj— i<f>{t) J x \ we obtain a 

Hamiltonian equivalent to (^|) 

H(t) = £JJ + [a ~ j>(t))J z + 2fl(t)j x . 

IV. SUMMARY 

In the present paper we have proposed and analyzed 
a robust adiabatic scheme for generating symmetric en- 
tangled states of a many-particle system. A constant 
nonlinear interaction in the collective spin projection J 2 
combined with a time-dependent linear interaction in J z 
creates a web of level crossings between the collective 
multiparticle states. Application of pulsed fields along 
the x component of the collective spin at the times of 
appropriate level crossings induces single or multi-step 
transitions between the states. In this way a controlled 
adiabatic navigation in the (N+ l)-dimensional subspace 
of the symmetric collective states is possible allowing, for 
example, generation of A^-particle W and GHZ states. 
The suggested method is robust against parameter vari- 
ations and does not request an exact knowledge of the 
particle number, but only demands proper timing of the 
pulsed field. It should be noted that the nonlinear inter- 
action J 2 guarantees the separation of the level crossings 
and is hence necessary for the navigation in the Hilbert 
space. 

An initial product state is connected to symmetric en- 
tangled states via various pathways involving avoided 
crossings induced by single and/or multiparticle inter- 
actions. If the exact number of particles is known, it 
is in general possible to find a pathway that runs only 
through single-excitation crossings; then only a moder- 
ately strong field is needed to induce adiabatic passage. If 
N is not known, a pathway can be chosen that involves 
multi-photon avoided crossings. Our numerical studies 
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indicate that in this case the required field scales nearly 
linearly with N. 

We have given two explicit examples for the applica- 
tion of our scheme, ion traps and coherently coupled two- 
component Bose condensates. 
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